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Interplay between changes of energy levels and tunneling amplitudes caused by localized electrons 
on-site Coulomb interaction depending on non-equilibrium electron filling numbers is analyzed. 
Specific features of local tunneling conductivity spectra for different positions of localized states 
energy relative to the Fermi level have been investigated by means of self-consistent mean field 
approximation in the presence of non-equilibrium effects. The conditions when modifications of 
tunneling transfer amplitude due to changes of electron filling numbers in the presence of on-site 
Coulomb interaction should be taken into account in tunneling conductivity spectra have been 
revealed. 
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I. INTRODUCTION 

Recently great attention is paid to both experimen- 
tal and theoretical investigations of kinetic processes 
in nanostructures, in connection with possibility of ul- 
tra small size electronic devices fabrication. Compari- 
son between experimental results and theoretical mod- 
els provide information about tunneling mechanisms in 
nanoscale junctions due to the presence of localized states 
[l| . Localized states formed by impurities at surface and 
interfaces of semiconductors lead to strong modification 
of local electronic structure and influence the behavior of 
tunneling characteristics due to on-site Coulomb inter- 
action of the sample and tip conduction electrons with 
charged localized states 0t3- Most of the experiments 
are carried out with the help of scanning tunneling mi- 
croscopy/spectroscopy technique ^2., ^ii] and theoretical 
calculations usually deals with Green's functions formal- 
ism [l|, or ab-initio calculaions • 

It was predicted theoretically and proved experimen- 
tally that on-site Coulomb repulsion of localized elec- 
trons leads to the changing of energy values considerably 
even for deep impurity levels [l|, Hf. It was found that 
taking into account Coulomb interaction of conduction 
electrons with non-equilibrium localized charges result 
in non-trivial behavior of tunneling characteristics in the 
case of STM metallic tip positioned above the impurity 
atom 0, [m HI] and also leads to formation of peculiari- 
ties in local tunneling conductivity analyzed apart from 
the impurity [1, . But in all the previous calculation 
tunneling transfer rates considered as a constant values 
and influence of the impurity's non-equilibrium electron 
flUing numbers on tunneling transfer amplitudes due to 
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Coulomb correlations was neglected |6H8l |12||. However, 
this effect is a problem of great interest because influence 
rate depends on the impurity type and solving of this 
problem gives an opportunity to initialize impurity type 
with the help of STM technique. Another question of 
great interest is what effect dominates in tunneling con- 
tact: influence of Coulomb interaction on-site or on hop- 
ping matrix elements. Correct investigation of this prob- 
lem shows the way for not only qualitative but also for 
quantitative theoretical analysis of experimentally mea- 
sured local tunneling conductivity, but it is necessary to 
choose the correct starting approximation for each type 
of impurity atom depending on the localized state en- 
ergy level position relative to the Fermi level. Previous 
works demonstrate that for localized states with energy 
deep below the Fermi level Hubbard-I model jT!] suites 
well and for localized states with energy in the vicinity of 
Fermi level mean field mixed valence approximation |14| 
can be applied. 

So in this work we present detailed investigation of 
local tunneling conductivity in vicinity of an impurity 
paying special attention to the tunneling transfer am- 
plitude dependence from non-equilibrium electron filling 
numbers due to Coulomb correlations. On-site Coulomb 
interaction will be also taken into account in our model. 
Development of theoretical description suggested in [§| 
and based only on on-site Coulomb interaction should 
clarify the role of tunneling transfer amplitudes modifi- 
cation due to the infiuence of impurity's non-equilibrium 
electron filling numbers in formation of peculiar features 
observed experimentally in local tunneling conductivity. 
We performed calculations for two extreme cases when 
tunneling takes place through the impurity state with 
deep energy level and shallow energy level. We applied 
modified mean field self-consistent approach for the im- 
purity state with shallow energy level and in the case 
of impurity with deep energy level modified Hubbard-I 
model was used. We also compared our results with the 
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calculations based on non-modified models based only on 
on-site Coulomb interaction. Both approaches were ana- 
lyzed with the use of the Keldysh fomalism [l^ . 



II. THE SUGGESTED MODEL AND MAIN 
RESULTS 



We shall analyze tunneling between the semiconductor 
surface and metallic STM tip in the presence of Coulomb 
interaction for impurities with deep and shallow energy 
levels in comparison with Fermi level. Special atten- 
tion will be payed to the impurity atom non-equilibrium 
electron filling numbers influence on the values of tun- 
neling transfer amplitudes. The standard model system 
semiconductor-impurity atom-metallic tip with constant 
tunneling transfer rates can be described by the Hamil- 
tonian: H: 



H — Hq + H, 
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Indices k and p label the states in the left (semicon- 
ductor) and right (tip) lead, respectively. The index d 
indicates that impurity electron operator is involved. U 
is the on-site Coulomb repulsion of localized electrons, 
iT-da = c^a'^da, Cda dcstroys impurity electron with spin 
a. Ho describes conduction electrons of continuous spec- 
trum states in the leads of tunneling contact. Htun de- 
scribes tunneling transitions between the tunneling con- 
tact leads through the impurity state. Himp corresponds 
to the electrons in the localized state formed by the im- 
purity atom. 

We shall start from the case of impurity state with en- 
ergy level situated deep below the Fermi level. In the 
presence of Coulomb interaction tunneling through deep 
impurity energy level can be described with the help of 
Hubbard-I model Ts"] (Hamiltonian of this model has the 
form [T]). Taking into account Coulomb interaction in 
this model leads to formation of two well separated im- 
purity energy levels Ed and Sd + U instead of one initial 
level Ed. It is reasonable to use approximation in which 
the strongest interaction of the considered model - the 
on-site Coulomb repulsion U - is included in G°f'^(a;). 
So we can write down expression for G'^^'^{uj) supposing 
tunneling transfer amplitudes to be independent from lo- 
calized state non-equilibrium filling numbers: 



UJ - Ed - S(w) 
nd-gUjuj - Ed) 
- Ed - (l - nd-cr)U 



(2) 



Using system of Dyson equations one can get expres- 
sion for impurity retarded Green's function. 



GR<7 
dd 



- Ed - S(w) - i(7fe<j -I- 7pd) 



(3) 



Where EptU^G^p? = Ipd and EAdl^Gl^ = l^d- 
Expression for local tunneling conductivity in the case of 
constant transfer amplitudes was found in [T7| and has 
the form: 



^( ) ^ Ji^j„,G'«f (w) 
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If one wants now to analyze modification of tunneling 
transfer amplitudes due to the dependence from impu- 
rity's non-equilibrium electron filling numbers, it is nec- 
essary to rewrite part Htun of the initial Hamiltonian [T] 
which corresponds to tunneling transitions between the 
tunneling contact leads through the impurity state. Due 
to Coulomb repulsion energy level position depends on 
electron filling numbers. Impurity energy increases with 
the increasing of filling numbers and it leads to the de- 
creasing of tunneling amplitudes between tunneling con- 
tact leads through the impurity energy level. So higher 
energy level corresponds to the smaller value of tunneling 
transfer amplitude and tunneling amplitude now depends 
on electron presence or absence on the energy level. Ini- 
tial Hamiltonian [1] can be rewritten in the following way: 



k,d(T p,da 

where expressions for t'^JJ' and t^^*^ have the form: 



tl" - < ^d-a >)+T2< Ud-a >) 

(tl(l- < Ud-a- >) +t2 < rid-a >) 
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pd 
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When impurity energy level is free from electrons tun- 
neling transitions from semiconductor continuous spec- 
trum states to STM tip through impurity atom are de- 
scribed by transfer amplitudes Ti and ti. If impurity 
energy level is occupied by electron with spin tunneling 
transitions between continuous spectrum states through 
impurity state can be described by tunneling transfer 
amplitudes T2 and ^2- Due to this assumption initial 
Hubbard-/ model [T3| can be modified and expression 
for retarded Green's function will have the form: 
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In this case relaxation rates can be written in the 
following way -fl{^f = Ektl^^rnGl^ and j^// = 

^ptp^J ImG'^^. Expression for local tunneling conduc- 
tivity than will have the form: 



dl -y-ff^-ff 



(8) 



Calculating expression for ImG^J {uj) requires solu- 
tion of self-consistent system of equations. Two of them 
are equations for ImG^^'^ iuj) , obtained from expression 
[71 and three equations determine impurity atom non- 
equilibrium electron filling numbers: 
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Ikd 



-Ipd 



(9) 



where n°(a;) and np{uj) are equilibrium filling numbers 
in the tunneling contact leads. 




FIG. 1: Local tunneling conductivity as a function of ap- 
plied bias voltage calculated for the impurity atom with deep 
energy level. Gray line corresponds to the case when tunnel- 
ing transfer amplitudes depend on impurity filling numbers 
ti = 0,20, t2 = 0,29, Ti = 0,30, T2 = 0,35. Black line 
corresponds to the case when tunneling transfer amplitudes 
have constant values and only on-site Coulomb interaction is 
taken into account ti^d = 0, 35, tpd = 0, 29. For all the figures 
parameters values Sd ~ —0, 80, U = 3 are the same. 

Fig. I shows tunneling conductivity as a function of ap- 
plied bias voltage for impurity with energy level deep be- 
low Fermi level. It means that relation between energy 
level value and energy level width due to the tunneling 
processes is large than unit. Tunneling conductivity cal- 
culated with the use of non-modified model Hubbard-I 



(Hamiltonian of this model has the form [T]) is shown 
by the black line and results obtained with the use of 
modified model (expression [5] for Htun) where depen- 
dence of tunneling transfer amplitudes from impurity 
non-equilibrium electron filling numbers is taken into ac- 
count are depicted by the gray line. Tunneling conductiv- 
ity calculated above the impurity demonstrates resonant 
peaks when applied bias voltage is equal to the initial 
position of impurity energy level (w = Sd) and when ap- 
plied bias voltage is equal to the initial position of impu- 
rity energy level shifted on the value of Coulomb poten- 
tial {uj — ed + U) both for taking (Fig. I gray line) and 
not taking into account dependence of tunneling transfer 
amplitudes from impurity non-equilibrium electron filling 
numbers (Fig. 1 black line). 

Modification of tunneling transfer amplitudes by non- 
equilibrium electron filling numbers results strongly in 
significant decreasing of resonant peaks amplitudes in 
comparison with the calculations performed with the 
use of Hubbard-I model for constant values of tunnel- 
ing transfer amplitudes (Fig. 1 black line). Resonant 
peculiarities also demonstrate considerable width chang- 
ing. When tunneling transfer amplitudes modification 
is taken into account resonant peaks width increases in 
comparison with calculations performed with the use of 
non-modified Hubbard-/ model. 

So one can clearly see that Coulomb interaction influ- 
ence on hopping matrix elements in addition to on-site 
Coulomb interaction strongly modify local tunneling con- 
ductivity measured above impurity with energy deep be- 
low Fermi level. Taking into account both effects leads 
to significant amplitude and width corrections of pecu- 
liarities in modified conductivity. 

Now let us analyze tunneling through the impurity 
state with energy level in the vicinity of Fermi level. It 
means that relation between energy level value and en- 
ergy level width due to the tunneling processes is about 
unit. In this case Coulomb interaction effects accompa- 
nied by the assumption that non-equilibrium electron fill- 
ing numbers modify tunneling transfer amplitudes can be 
described with the use of modified mean-field mixed va- 
lence approximation [3, consequently Hamilto- 
nian depends on average electron filling numbers. So one 
has to rewrite again part Htun of the initial Hamiltonian 
[l]in the form[5l More over impurity energy level position 
depends on Coulomb interaction of the non-equilibrium 
electron density. A new impurity energy level position 
is a result of initial energy level position shift on the 
value U- < Ud-a >■ Non-equilibrium impurity atom fill- 
ing numbers Ud-a must satisfy self-consistency condition 
[5] where ^"^^J and ^^^J are determined by expressions 

[6] and G^^{uj) in the mean field approximation has the 
following form: 



dd 



1 



(10) 



After determining new position of impurity energy 
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level we can calculate expressions for ImG^^ (uj) from 
1101 and finally evaluate local tunneling conductivity for 
shallow impurity taking into account Coulomb interac- 
tion effects and tunneling transfer amplitudes modifica- 
tion influenced by the impurity non-equilibrium electron 
filling numbers with the use of expression |8l 

0,5- 



0,4- 




(0 ,eV 

FIG. 2: Local tunneling conductivity as a function of applied 
bias voltage calculated for the impurity atom with shallow 
energy level. Gray line corresponds to the case when tunnel- 
ing transfer amplitudes depend on impurity filling numbers 
ti = 0,40, t2 = 0,49, Ti = 0,60, T2 = 0,68. Black line 
corresponds to the case when tunneling transfer amplitudes 
have constant values and only on-site Coulomb interaction is 
taken into account tkd = 0, 68, tpd = 0, 49. For all the figures 
parameters values Sd ~ —0, 80, U = 3 are the same. 

Calculation results are presented on Fig. 2. Tunneling 
conductivity as a function of applied bias voltage demon- 
strates that taking into account modification of tunnel- 
ing transfer amplitudes due to the influence of impurity's 
non-equilibrium electron filling numbers (Fig. 2 gray line) 
in addition to Coulomb interaction leads to important 
changes in comparison with the case of constant tun- 
neling rates and on-site Coulomb interaction(Fig.2 black 
line) . Decreasing of resonant peak width and a negligible 
decreasing of the tunneling conductivity amplitude takes 
place for modified self-consistent mean field approxima- 
tion due to the influence of impurity's non-equilibrium 
electron filling numbers on tunneling transfer amplitudes 
and consequently on relaxation rates. Describing tun- 
neling through the impurity state with energy level in 
the vicinity of Fermi level it is necessary to mention 
that width of peculiarities in tunneling conductivity is 
determined by the value of relaxation rates 7^!^''^ and 
Tpd^- Continuous energy shift of the resonant peak to 
the higher values of applied bias voltage exists mostly 



due to on-site Coulomb interaction of localized electrons 
{U < Ud^a »"f'^^^) and can be seen both for taking and 
not taking into account tunneling matrix elements modi- 
fication. It is clearly evident that resonant peak shift for 
modified tunneling amplitudes is smaller than in the case 
of constant tunneling amplitudes also because of relax- 
ation rates modification. 

III. CONCLUSION 

In this work we have analyzed the role of tunneling 
transfer amplitudes modification caused by the shift of 
energy level with changing of non-equilibrium electron 
filling numbers in the presence of on-cite Coulomb inter- 
action in formation of anomalous features in local tun- 
neling conductivity. The influence of this effect on local 
tunneling conductivity have been analyzed for different 
location of localized states energy relative to Fermi level. 
We have studied two extreme cases when tunneling takes 
place through the impurity state with energy level posi- 
tion deep below Fermi level and through impurity state 
with energy level in the vicinity of Fermi level by means 
of modified Hubbard-I model and modified self-consistent 
mixed valence mean field approximation correspondingly. 

For impurity state with energy level deep below the 
Fermi level modification of tunneling transfer amplitudes 
in addition to on-site Coulomb interaction results in sig- 
nificant changing of resonant peaks amplitudes and width 
in comparison with non-modified Hubbard-I model. Res- 
onant peculiarities demonstrate strong amplitudes de- 
creasing. We also observed considerable increasing of 
resonant peculiarities width. 

Impurity state with energy in the vicinity of Fermi level 
also demonstrates significant changes in local tunneling 
conductivity in comparison with the case of constant tun- 
neling rates due to modification of tunneling transfer am- 
plitudes in addition to on-site Coulomb interaction. We 
observed width changing of resonant peak, slight decreas- 
ing of the tunneling conductivity amplitude and smaller 
resonant peak shift to higher value of applied bias volt- 
age. 

Obtained results give us possibility to conclude that 
on-site Coulomb interaction strongly influence local tun- 
neling conductivity measured in the vicinity of impurity 
atom and Coulomb interaction influence on hopping ma- 
trix elements leads to corrections of peculiarities shape, 
amplitude and width in modified conductivity. 

This work was supported by RFBR grants and by the 
National Grants for technical regulation and metrology. 
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